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ABSTRACT 

We calculated the Knight shift and spin-lattice relaxation rates of Pu-115 compounds assuming 
O d-wave superconductivity in the presence of strong impurity scattering. We discuss the implica- 
^ ■ tions for recent measurements of the spin- lattice relaxation rate in the Pu-115 compound PuRhGas 
^ ■ by Sakai and coworkers [J. Phys. Soc. Jpn. 74, 1710 (2005)] and present a prediction for the 
^ ! corresponding Knight shift. In addition, we noticed a significant round-off of the spin-lattice re- 
\ laxation rate 1/Ti just above the superconducting transition temperature that is not observed 
£h ■ in the sister compound PuCoGas. It appears that in PuRhGas superconductivity is mediated 
^ij by spin fluctuations, too. This provides additional support to the scenario of superconducting 
£H [ pairing mediated by spin fluctuations in the Pu-115 compounds similar to the Ce-115 compounds 
O \ and the high-temperature copper-oxide superconductors. 

T 1 ■ INTRODUCTION 

> ■ 

The discovery of superconductivity in plutonium based systems such as PuCoGa 5 1] and 
PuRhGa 5 !2j has stimulated the study of unconventional superconductivity and the pairing sym- 
t-h ■ metry and mechanism in these materials. The symmetry of an unconventional superconductor is 
^ . reduced compared to the symmetry of its normal state, thus resulting in many novel properties 
-j-^. \ of the quasiparticle excitation spectrum. It is believed that the superconducting action in Pu- 
0$ | 115 [PuMGas with M=Co and Rh] derives itself from the unique character of the 5f electrons of 
plutonium [3|. The tetragonal crystal structure of PuMGa5 is isostructural to that of the Ce-115 
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^ . series [CeMIn, 

The purpose of this study is to shed light on the superconducting pairing symmetry and possible 
^ ' pairing mechanism in the Pu-115 compounds. Very recently, Curro and coworkers proposed, 
based on their measurements of the Knight shift and spin-lattice relaxation rates, that the Pu-115 
k~j com P oun ds are bridging the superconducting and normal-state properties of the heavy-fermion Ce- 
rN ' 115 and high-temperature copper-oxide superconductors. Therefore providing a means for tuning 
& ' the interaction strength of antiferromagnetic spin fluctuations to intermediate values between 
both extreme limits [5J]. 

The experimental techniques of nuclear magnetic resonance (NMR) and nuclear quadrupolar 
resonance (NQR) have been used successfully in the past to distinguish between the spin states 
of Cooper pairs (spin singlet vs. spin triplet pairing) and provide indirect information on the 
symmetry of the gap function - fully gapped vs. nodal lines or nodal points in the gap function 
on the Fermi surface. Both techniques probe directly the quasiparticle density of states and reveal 
indirect information about the pairing symmetry. 

The standard explanation of power vs. exponential laws in the low-temperature behavior of 
thermodynamic and transport properties, for example, the spin-lattice relaxation rate 1/Ti, comes 
from the difference of nodal and fully gapped excitation spectra in the superconducting state. 



In clean nodal superconductors 1/T% exhibits a nearly T 3 behavior far below the superconduct- 
ing transition temperature T c , while it is exponential for gapped superconductors. On the other 
side, deviations from this behavior, like the T-linear temperature dependence of 1/T\ at low tem- 
peratures, are explained by impurity effects in an unconventional superconductor (SC) with lines 
of nodes on the Fermi surface. Very recently, Sakai and coworkers p reported such a result for 
the spin-lattice relaxation rate of PuRhGas in the superconducting phase. This behavior closely 
resembles the spin-lattice relaxation rates measured by Curro and coworkers for PuCoGa^, which 
belongs to the same family of Pu-115 compounds, albeit with a T c nearly twice as high |^|. 

Here we give a detailed theoretical description of the spin-lattice relaxation rate and predict 
what should be observed for the Knight shift if measured on the same sample. Our self-consistent 
treatment of impurity scattering in the superconducting state goes beyond the two-fluid approach 
used by Sakai et al. |.6;], which was used to explain the large residual density of states in PuRhGa 5 . 
Simultaneous measurements of spin-lattice relaxation rate and Knight shift will place stringent 
constraints on the symmetry and magnitude of the superconducting gap function, as well as on 
the concentration and scattering strength of impurities at low temperatures. 



THEORY 

The effect of impurity scattering is included within the self-consistent T-matrix approximation 
7, 8|, which is the standard formulation for pointlike defects in a superconducting dilute alloy 
Si llO, III |- For the case of particle-hole symmetry of the quasiparticle excitation spectrum the 
Nambu component T3 of the T matrix vanishes, and for a d-wave order parameter (OP) with 
isotropic scattering T\ = (also without loss of generality we can choose T2 = by general U(l) 
gauge symmetry), where Tj is the ith component of the 2x2 Nambu matrix expanded in Pauli 
matrices. Then we need to calculate only T (u). The impurity self-energy is given by S = TT , 
where T = rii/nNQ. Here N Q is the normal density of states (DOS) at the Fermi surface (FS), 
m is the impurity concentration; T (u n ) = [e2 !_°^ ( "J n)] , where g (u n ) = ^ Efc a2 + ffi A 2 (fc) - The 
impurity renormalized Matsubara frequency is defined by u n = u) n + E , with u n = nT{2n + 1), 
and the scattering strength parameter c is related to the s-wave phase shift S by c = cot(<5 ). 
Using this self-energy S the following gap equation is solved self-consistently, 



A(0) = -N g(<f>) [ _^y(0-0')T]T f C de^ 
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where V(<f) — <f)') is the angular parametrization of the pairing interaction, and u c is a typical cutoff 
energy. We assume the canonical d-wave gap function of the form A (A;) = A (cosfc x — cosk y ) 
or A(0) = Aocos(20) for a cylindrical Fermi surface. The pairing potential V(<f) — ) induces 
a gap with d-wave symmetry. Although its microscopic origin is not the issue of this paper, we 
believe it originates from antiferromagnetic (AFM) spin fluctuations. The static limit of the spin 
susceptibility of the AFM fluctuations, xi^ 00 = 0) ~ (q-Q^+g- 2 > ^ s parameterized near the AFM 
wave vector Q as [l2| 

n0-0') = K, W (0 _ , ± ^ /2)2+fc2 , ( 2) 

where the parameter b is inverse proportional to the AFM correlation length £, normalized by the 
cylindrical FS (£ ~ arr/b; a is the lattice parameter). For all calculations in this paper, we chose 
b = 0.5 which is not a sensitive parameter for our results unless £ is very large (6 < 0.1) [~L2j] . i.e., 
within the range of 0.1 < b < 1 our results show little variations and are qualitatively the same. 



With the gap function A(0) and To (a;) obtained from Eq. (JTJ) (T (uj) is analytically continued 
from T Ujn) byPade approximant method), we calculate the l/Ti nuclear spin-lattice relaxation 
rate 0,0, 0,0 
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and the superconducting spin susceptibility xs 
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where is the Fermi-Dirac function, the impurity renormalized quasiparticle energy uj = 

uj + Hq(u), and (...)$ means the angular average over the FS. The first term in the bracket of 
Eq. (JHJ) is N 2 (u). The second term vanishes in our calculations because of the symmetry of the 
OP. To calculate \jT\T using Eq. (J3J), or xs using Eq. (jlj), we need the full temperature dependent 
gap function A(0, T) and T c . Our gap equation Eq. (JTJ) is the BCS gap equation, therefore it gives 
the BCS temperature behavior for A(0, T) and A = 2.14 k B T c for the standard weak-coupling 
d-wave SC. In order to account for strong-coupling effects we use the phenomenological formula 



A(0, T) = A(0, T = 0) S(T) with H(T) = tanh(/3y T c /T — 1), and parameters (3 and A /T c . Then 
we only need to calculate A(0, 0) at zero temperature. The temperature dependence of So(co>,T) 
(= T T (u, T)) is similarly extrapolated: T (u, T) = T (u, T = 0) H(T) + T norma i(l — E(T)), where 
T nor mai = P/(c 2 + 1) is the normal state T . In our numerical calculations we chose (3 = 1.74, 
because our final results are not very sensitive with respect to this parameter, while the ratio 
A /fcsT c is an important parameter to simulate strong-coupling effects. The larger this ratio is, 
the more important are strong-coupling effects. 



RESULTS AND DISCUSSIONS 

In figures ^ and El the spin-lattice relaxation rate of PuRhGas by Sakai et al. is shown, 
where l/Ti is normalized to l/Ti = 10 at T = T c for ease of comparison. The insets show 
the corresponding normalized quasiparticle DOS for varying scattering rates T. With our earlier 
described choice of parameters the impurity scattering rate T/Aq = 0.032 is enough to completely 
fill the low energy gap with impurity states and N(uj = 0) reaches more than 25% of the normal- 
state DOS iVo- 

In Fig. n we obtain a better fit to the experimental data (symbols) assuming a slightly lower 
superconducting transition temperature T c = 7.6 K than their reported value of T c = 8.5 K. 
This could indicate the presence of a pseudogap similar to the high-temperature superconductor 
YBa2Cu307_,5, where l/Ti is suppressed just above T c . 

For the temperature dependence of the gap, we chose the parameters f3 = 1.74 and 2A = 
5k B T c for the d-wave gap to account for the strong-coupling effects of superconductivity as ex- 
plained before. As expected from the DOS results, due to the impurity induced residual states, 
l/Ti displays the linear-T dependence at low temperatures and the region of T-linear behavior 
increases with impurity concentration. For a higher value of r/A = 0.064, this T-linear region 
extends up to ~ 0.35 T c . At temperatures near T c the coherence peak is almost invisible because 
of the sign-changing gap function, i.e., vanishing of the second term in Eq. (jSJ). Below T c it shows 
a nearly T 3 behavior due to the lines of nodes in the gap until it goes through a gradual crossover 




FIG. 1: The NQR spin-lattice relaxation rate plotted versus temperature normalized by T c . Calculations 
are for 2Aq = hk^T c and three values of the impurity scattering rate T for unitary scattering. Inset: 
The normalized quasiparticle DOS for corresponding values of T/Aq = 0,0.032,0.064. 




FIG. 2: The NQR spin-lattice relaxation rate plotted versus temperature normalized by T c . Calculations 
are for 2Ao = 8&bT c and three values of the impurity scattering rate T for unitary scattering. Inset: 
The normalized quasiparticle DOS for corresponding values of T/Aq = 0,0.032,0.064. 
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FIG. 3: The calculated spin susceptibility xs of a d-wave SC normalized by its normal state value xn 
for gap values 2A = 5/cbT c (solid lines) and 8&bT c (dotted lines), and impurity scattering T/A = 
and 0.032. 

region and finally to the T-linear region. The comparison with the experimental data by Sakai et 
al. on PuRhGas is in good agreement with unitary scattering, a gap value 2A = 5 /c_bT c , and 
a scattering rate close to r/A = 0.032. Based on this value, we estimate the superconducting 
transition temperature of the pristine sample to be T c0 = T c + ~T « 8.1 K or 9.0 K, depending 
on the value of T c = 7.6 K or 8.5 K. 

In Fig. |21 the normalized 1/Ti is plotted for an enhanced strong-coupling d-wave gap value 
2A = 8 k B T c , as was recently found for PuCoGa 5 0|. Due to the larger gap value, the calculated 
1/Ti is always less than the measured spin-lattice relaxation rate. Hence we find a poorer fit to 
the experimental data for this choice of the strong-coupling gap. 

Fig. El shows the prediction for the spin susceptibility, xs, or its corresponding NMR Knight 
shift, K = Kq + A\si where K and A are constants for most materials, xs is calculated for the 
same d-wave gap values as was used for the spin-lattice relaxation rates in figures ^ and El Again 
a modest impurity scattering rate of T/Aq = 0.032 results in a large residual susceptibility at 
zero temperature, equivalent to roughly 25% of the normal state DOS or spin susceptibility xn- 
The quantitative difference in the spin susceptibility between gap values 2A = 5 ksT c and 8 fcsT c 
should be easily discernible in measurements of the Knight shift. 

CONCLUSIONS 

The NQR spin-lattice relaxation rate 1/Ti in PuRhGas is consistent with a strong-coupling 
d-wave gap function and unitary impurity scattering similar to the observed behavior in its sister 
compound PuCoGas. Inspite of many similarities between PuCoGas and PuRhGas, we also find 



marked differences: (1) The maximum superconducting gap value of PuRhGas is smaller than for 
PuCoGa 5 , i.e., it is 2A /A;bT c = 5 for PuRhGa 5 versus 8 for PuCoGa 5 . This suggests that the 
fluctuations of the pairing bosons are weaker for PuRhGas, possibly hinting at a progressive trend 
for the strength of the spin-fluctuations in this class of materials. (2) Although the PuRhGas sam- 
ple was of similar age as PuCoGa 5 when measured, it had a three times larger relative scattering 
rate r/A = 0.032 compared to PuCoGa 5 with r/A = 0.01, which could be due to variations 
of the isotope mix of plutonium between both samples. (3) 1/Ti exhibits a rounded behavior 
between T = 7.6 K and 9 K resembling the pseudogap phenomenon in YBa2Cus07_5, while no 
such rounding is observed for PuCoGas. This certainly needs clarifications by further studies of 
1/Ti in the normal state, which will be reported in a separate paper. 
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